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Predicting the ground states for surface adsorption is a challenging problem because the number 
of degrees of freedom involved in the process is very high. Most of the studies deal with some 
specific arrangements of adsorbates on a given surface, but very few of them actually attempt to 
find the ground states for different adatom coverage. In this work, we show the effectiveness of cluster 
expansion method to predict the “ground states” resulting from chemisorption of oxygen and fluorine 
atom on the surface of monolayer black phosphorus or phosphorene. For device applications, we 
find that in addition to bandgap tuning, controlled chemisorption can change the unique anisotropic 
carrier effective mass for both the electrons and holes and even rotate them by 90°, which can be 
useful for exploring unusual quantum Hall effect and novel electronic devices based on phosphorene. 

PACS numbers: 68.43.-h, 73.22.-f, 73.20.At 


I. INTRODUCTION 

Black phosphorus has emerged as a very promising 2D 
material, on account of its direct band gap, high car¬ 
rier mobility, thickness dependent and highly anisotropic 
electronic properties.Moreover, few layers of black 
phosphorus has also been successfully exfoliated from its 
bulk counterpart and is available in free standing form.^“^ 
As a shortcoming, due to the presence of a lone pair 
of electrons, phosphorus is highly reactive at the sur¬ 
face, and this can lead to the degradation of phospho¬ 
rene (for example, via photo oxidation^”^) under ambient 
conditions, hindering its potential applications. A 
possible solution to this problem may be chemical pas¬ 
sivation of phosphorene with adsorption of monovalent 
and divalent atoms such as fluorine and oxygen. Apart 
from controlled passivation of phosphorene, chemisorp¬ 
tion may also lead to new stable structures with inter¬ 
esting electronic properties as in the case of graphene, 
silicene and 

As observed by Ziletti et al.^^ based on first principles 
calculations, chemisorption of a single oxygen atom at the 
phosphorene surface is an exothermic process. Interest¬ 
ingly, dangling or surface location (O bonded to only one 
P atom) is energetically more favorable than the bridge or 
interstitial position (O bonded to two P atoms).Such 
a defect, where O is dangling, causes minimal lattice dis¬ 
tortion and it is also electrically neutral with no midgap 
states.Not only the bridge position is energetically 
less favorable because of large lattice distortion, but also 
it requires 0.7 eV of activation energy for one O atom to 
hop from a dangling to a bridge location.In case of high 
O concentration (upto 250%), various planar and tubular 
forms of oxides have been proposed with the former being 
energetically more favorable.For the case of monova¬ 
lent atom adsorption, first principles calculations predict 
better stability for fluorinated phosphorene as compared 
to hydrogenated phosphorene, which is found to be dy¬ 


namically unstable. 

One important but unexplored problem is to find the 
ground states among the many possible arrangements of 
adsorbate on phosphorene for a given adatom concen¬ 
tration. Since it involves sampling of a 2^ dimensional 
configuration space {N is the number of lattice sites oc¬ 
cupied by P atoms, each of which can either be coupled 
to an adatom or remain vacant), a ground state search 
by exhausting all the possible combinations by means for 
first principles calculations is not feasible. As a feasible 
alternative, we choose cluster expansion (CE) technique 
for searching the ground state at various different concen¬ 
tration in between pristine and fully oxidized/fluorinated 
phosphorene (PO/PF). Unlike mean field approximation, 
the CE technique is well known for giving a complete mi¬ 
croscopic description of the system. 

In this article, we predict several new “ground states” 
of fluorinated and oxidized phosphorene and also study 
the effect of adsorption on electronic properties with 
varying concentration of the adsorbate. Note that, we 
consider chemisorption only at the dangling position lo¬ 
cated at the surface and these are possibly not the univer¬ 
sal ground states, in particular at higher concentration of 
divalent oxygen, which can also be present at the inter¬ 
stitial bridge position. However, we argue later that such 
metastable states of oxidized phosphorene, predicted in 
this paper, can be stabilized kinetically and thus, pos¬ 
sibly be synthesized under suitable experimental condi¬ 
tion. Since fluorine is monovalent, it can only be ad¬ 
sorbed in the dangling position and we believe that the 
universal ground states are correctly predicted in this 
work. In terms of electronic properties, we find that 
other than bandgap modification, systematic function¬ 
alization of phosphorene can tune the anisotropy of both 
the electron and hole effective masses in a controlled way; 
unlike the case of external strain, which can tune only 
the electron effective mass.^’^® We also check the thermal 
stability of all the “ground state” structures by doing a 
molecular dynamics (MD) simulation. This article is or- 
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ganized as follows. In Sec. II, we discuss the details of 
our computational approach. This is followed by a thor¬ 
ough discussion of fluorination of monolayer phosphorene 
in Sec. Ill, and a discussion of oxidation of monolayer 
phosphorene in Sec. IV. Thermal stability of discovered 
“ground states” are discussed in Sec. V and finally we 
summarize our results in Sec. VI. 


II. COMPUTATIONAL DETAILS 

The ground state search at various different adatom 
concentrations is performed by using cluster expansion 
(CE) technique,^®^^^ as implemented in alloy-theoretic 
automated toolkit (ATAT).^°^^^ As the name suggests, 
the code can predict intermediate ground states of bi¬ 
nary alloys, made by mixing different fractions of two 
pure components A and V. Surface adsorption can be 
treated using this method via a mapping onto the pseudo¬ 
alloy problem by first identifying adatom locations (in the 
phosphorene unit cell), which can either be occupied by 
the adatom (component A) or remain vacant (component 
V). Thus, pristine and fully fluorinated/oxidized phos¬ 
phorene corresponds to 100% V - 0% A and vice versa, 
respectively. CE is a technique to parameterize the total 
energy for any given configuration of A and V for an ar¬ 
bitrary composition, expressed as (0 < a: < 1), 

thus avoiding the cumbersome task of computing total 
energy of 2^ (N is of the order of lattice sites) structures 
from first principles calculations. As per the CE formal¬ 
ism, any given configuration of A and V can be expressed 
by the vector cr = {ai,U 2 , where ai = -1-1 or —1 

depending on whether site is occupied by A (F/0) or 
V (vacant). Then, energy for a given configuration cr is 
expressed by an “Ising-like” Hamiltonian, 

Eicr) ='^maJa 

a Xiep / 

where a represents a cluster, which is a collection of sites 
j, for example pairs, triplets, quadruplets etc. The multi¬ 
plicities m^ is the number of clusters of type a and Jq, is 
the corresponding effective cluster interaction (ECI). The 
sum is carried out over all clusters a, which are not equiv¬ 
alent by space group symmetry and the “spin” product is 
averaged over all the clusters /3, those are equivalent to a- 
type clusters by symmetry. The ECIs (J^) are obtained 
by fitting the energies (of selected structures, generally 
50-100 in number), obtained directly from ah initio cal¬ 
culations based on Density Functional Theory (DFT).^^ 
The quality of the fit is determined by calculating the 
cross-validation score, 

n 

= ( 2 ) 

where is the DFT calculated energy of a particu¬ 

lar structure i and Ei is the predicted energy of the same 


structure, calculated from a least-square fit to the ener¬ 
gies of the rest of the (n-1) structures (excluding E^’^'^). 

Once the energy is parameterized in terms of ECIs 
(Ja), ground state for a particular adatom concentration 
X is determined by comparing the formation energies, de¬ 
fined as 

Ef = E{A,Vi.,) - xEiAiVo) - (1 - x)EiAoVi) (3) 

where x denotes the adatom concentration and E{AiVq) 
and E{AqVi) are the energies of fully functionalized and 
pristine phosphorene, respectively. Relative stability of 
all the predicted ground states is compared in terms of 
binding energy, given by 

Eb = [Ep+a — {EfpEp + NaEa)] (4) 

where Na and Np denote the number of adatoms and 
phosphorus atoms in the unit cell, respectively. In 
Eq. (4), Ep+a, Ea and Ep are the energies of adatom 
adsorbed phosphorene monolayer, a single phosphorus 
atom in pristine phosphorene monolayer and a free atom 
of fluorine or oxygen (depending on whether we consider 
fluorination or oxidation), respectively. 

The first principles (DFT) calculations are per¬ 
formed using a plane-wave basis set and ultrasoft 
pseudo-potential, as implemented in Quantum Espresso 
package.Electron exchange-correlation is treated with 
a generalized gradient approximation (GGA), using 
the formulation developed by Perdew-Burke-Ernzerhof 
(PBE).^^ Kinetic energy cutoff for the plane-wave basis 
for wave functions is taken to be equal to 40 Ry. A well 
converged fc-point mesh, depending on the size of each 
unit cell (14 x 16 x 1 for the smallest unit cell of pristine 
phosphorene), is chosen for the Brillouin-zone integra¬ 
tions. The structures are relaxed (both cell-vectors and 
ionic positions) until the the forces on each atom (to¬ 
tal energy change due to ionic relaxation between two 
successive steps) is less than I0“® Ry/au (10“® Ry). A 
supercell with a vacuum of ^25 A perpendicular to the 
atomic plane is used to avoid the interaction with spu¬ 
rious replica images. Using these parameters, the in¬ 
plane and out-of-plane bond lengths of pristine mono- 
layer phosphorene are found to be 2.23 and 2.26 A, re¬ 
spectively, which are in very good agreement with the 
values reported in the literature®. Additionally we will 
see later that the binding energy of F/0 atoms is very 
high, which clearly shows that these atoms are chemically 
bonded to the P atoms, as a consequence we do not need 
to take into account van der Waals interactions for this 
problem. However, in case of multiple layers, it would be 
appropriate to include van der Waals interactions, which 
mainly affects the inter-layer spacing. 

III. FLUORINATION OF PHOSPHORENE 

Unlike silicene, which is grown on a suitable sub¬ 
strate,®® monolayer phosphorene is a 2d-material that 
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FIG. 1. (a) Formation energies (Ef) at different fluorine concentration, calculated from DFT (triangles) and fitted values (cross) 
from Eq. 2. Blue circles indicate the ground states, obtained for x = {0,0.25,0.5,0.75,0.833,0.875,1.0}, whose binding energy 
[calculated using Eq. (4)] and bandgap values are plotted in (b) and (c), respectively, (d) Unit cell (top and side view) of each 
of the ground states with increasing concentration of fluorine, and (e) the corresponding electronic band-structures, (f) The 
orbital resolved density of states (DOS) of the respective ground states, with the inset denoting the total DOS. Note that the 
solid and dashed lines correspond to phosphorus and fluorine, respectively. 


can be mechanically exfoliated^® and it is found to be 
stable in freestanding form, similar to graphene®^ and 
M 0 S 2 .®® Thus, our ground state search includes both sin¬ 
gle sided as well as double sided functionalization. How¬ 
ever, we find that the lowest energy structures are always 
those formed by double sided adsorption in case of fluo- 
rination, as well as oxidation. Fluorine is chosen as the 
monovalent adsorbate because hydrogenated structures 
of phosphorene are found to be unstable. 

Applying cluster expansion technique, we first deter¬ 
mine the ECI based on total energy calculation of 74 
structures using DFT method. A cross validation score 
of 20 meV reflects that the difference between fitted 
and actual (from DFT calculations) energy of the flu- 
orinated structures is negligible. Thus, based on the cal¬ 


culated ECI values, formation energies of 18652 fluori- 
nated phosphorene structures (up to 32 atoms per unit¬ 
cell) are predicted and the true and predicted ground 
states agree very well [see Fig. 1(a)]. We find six ground 
states of fluorinated phosphorene: P 8 F 2 , P 4 F 2 , PsFe, 
P 12 F 10 , P 16 F 14 and P 4 F 4 at specific fluorine concentra¬ 
tions of a; = 0.25,0.5,0.75,0.833,0.875 and 1.0, and the 
atomic arrangements, along with the unit cells are shown 
in Fig 1(d). The binding energy decreases with increas¬ 
ing fluorine concentration [see Fig 1(b)]. Note that, even 
at low F concentration the ground state configuration is 
constituted by two sided adsorption [see Fig 1(d)], which 
clearly indicates that one sided fluorination is energeti¬ 
cally less favorable. A clear pattern is also visible, where 
the adatom located on a phosphorus atom in the upper 
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FIG. 2. Polar plot of the hole and electron effective masses 
for various ground states of fluorinated phosphorene. Note 
that, has a two fold rotation symmetry for all the ground 
states, other than the fully fluorinated phosphorene, where 
the four fold rotation symmetry is evident. On the other 
hand, m^ maintains it’s two fold rotation symmetry, although 
becoming less anisotropic after complete fiuorination. 


plane of the puckered structure is accompanied by an¬ 
other adsorbed fluorine on the neighboring phosphorus 
atom in the lower plane. Fiuorination leads to significant 
distortions in the structural parameters of the underlying 
monolayer phosphorene, which in the pristine form has a 
bond-length of 2.2 A and bond-angles of 96° and 104°. 
The structural parameter data for pristine phosphorene 
are in good agreement with the values reported in the 
literature.^® On fiuorination, the bonds corresponding to 
the pair of phosphorus atoms to which the F atoms are 
adsorbed get elongated (bond-length ranging from 2.5 A 
to 3 A) while the bond-angles show small variations (up- 
to 5%). The remaining P-P bonds, to which fluorine is 
not attached show negligible distortions. 

Next we focus on electronic band structure, which 
changes significantly due to fluorine adsorption. Based 
on spin-resolved electronic structure calculations we find 
that fiuorination does not induce any magnetization in 
phosphorene. Figure 1(e) displays the electronic band 
structure of fluorinated phosphorene with increasing flu¬ 
orine concentration. The energy dispersion relation is 
plotted along the high symmetry directions (S-Y, Y-F, 
F-X, X-S and S-F), where the high-symmetry points in 
the Brillouin zone are as follows: r(0,0,0), X(0.5,0,0), 
Y(0,0.5,0) and S (0.5,0.5,0). The direct band gap of 
pristine phosphorene (0.9 eV at F point, within GGA) 
is initially found to decrease to 0.1 eV (Y point) at 
x=0.25. With further increase in fluorine concentration, 
the bandgap undergoes a transition from direct to indi¬ 
rect, accompanied by a monotonically increasing mag¬ 
nitude, reaching a value of 1.19 eV in fully fluorinated 
phosphorene (x=l). This variation of the bandgap as a 
function of fluorine concentration is shown in Fig. 1(c). 
Even though fiuorination significantly alters the elec¬ 
tronic structure of pristine phosphorene, both conduction 


band minima (CBM, corresponds to electron like charge 
carriers) and valence band maxima (VBM, corresponds 
to hole like charge carriers) states are still dominated by 
the orbitals belonging solely to the phosphorus atoms. 
Thus fluorine adatoms are not going to play any active 
role in carrier transport in black phosphorene. However, 
unlike pristine monolayer phosphorene where orbitals 
contributes most to the VBM states, it is the Px and Py 
orbitals which dominate in fluorinated phosphorene [see 
Fig. 1(f)]. In case of conduction band edge, s as well as 
all the p orbitals of phosphorus have significant weightage 
in pristine as well as fluorinated phosphorene. 

One of the exciting properties of phosphorene is it’s 
anisotropic electrical conductance and optical responses,^ 
which is attributed to the anisotropy of the effective 
mass of the charge carriers.As shown in Fig. 1(e), 
compared to the F-Y direction, both the valence and 
conduction band in the F-X direction are relatively flat 
in pristine phosphorene near the F point. Thus, both 
electron (m*) and hole (m)() effective mass is higher in 
the F-X direction (m*=1.2 mg and m)j=8.64 me, where 
mg « 9.11 X 10“®^ kg), than compared to their values 
in the F-Y direction (me=0.17 me and m)[=0.16 mg). 
The effective mass tensor is diagrammatically presented 
in Fig. 2, which approximately looks like the number 8 
for pristine phosphorene. It has been shown that by ap¬ 
plying strain, the electron effective mass tensor can be 
rotated by 90°, taking the shape of In contrast, 

the hole effective mass is not so sensitive to the external 
strain; the nature of the anisotropy remains unchanged, 
with significantly higher effective mass in F-X direction 
than compared to F-Y direction.^’®® From an application 
point of view, this phenomenon can be used to switch 
the easy direction of electronic and thermal transport in 
phosphorene based devices. 

Since fiuorination changes the location of VBM and 
CBM, accompanied by significant alteration of curvature 
of the respective bands [see Fig. 1(e)], we explore whether 
adatom adsorption can be an alternative way of tuning 
the anisotropic effective mass of phosphorene. As shown 
in Fig. 2, m)[ in the F-X direction decreases significantly 
from it’s value of 8.64 mg in pristine phosphorene, ly¬ 
ing in a range from 0.5 mg to 1.2 mg after fiuorination. 
On the other hand, m)[ in F-Y direction does not change 
significantly after fiuorination (ranging from 0.15 mg to 
0.25 mg, which is very close to it’s value of 0.16 mg, found 
in pristine phosphorene), except for high concentration of 
fluorine. As shown in Fig. 2, value of m^ in F-Y direction 
for x=0.875 and x=1.0 is 1.28 mg and 0.55 mg, respec¬ 
tively. Interestingly, only for x=0.875, the anisotropy 
of m)( is rotated by 90° and it’s value in F-Y direction 
is nearly two times higher than that of F-X. In case of 
fully fluorinated phosphorene, the effective mass plot has 
a butterfly like shape. Broadly speaking, fiuorination of 
phosphorene leads to significant decrease of m)) in the 
F-X direction, while it is weakly affected in the F-Y di¬ 
rection, except at high adsorbate concentration. 

In case of electrons, effective mass in the F-X direction 
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FIG. 3. (a) Formation energies (Ef) at different oxygen concentration, calculated from DFT (triangles) and fitted values (cross) 
from Eq. 2. Blue circles indicate the “ground states”, obtained for x = {0,0.167,0.25,0.33,0.5,0.75,0.917,1.0}, whose binding 
energy [calculated using Eq. (4)] and bandgap values are plotted in (b) and (c), respectively, (d) Unit cell (top and side view) 
of each of the “ground states” with increasing concentration of oxygen, and (e) the corresponding electronic band-structures, 
(f) The orbital resolved density of states (DOS) of the respective “ground states”, with the inset denoting the total DOS. Note 
that the solid and dashed lines correspond to phosphorus and oxygen, respectively. 


decreases due to fluorination, lying in a range from 0.1 mg 
to 0.3 mg, other than the case of fully fluorinated phos- 
phorene (1.1 mg), where it is nearly equal to m* of pris¬ 
tine phosphorene (1.2 mg). On the other hand, effective 
mass calculated along the F-Y direction increases mod¬ 
erately; after fluorination m* lies in the range from 0.25 
iTLg to 0.75 mg except for the case of x=0.25, for which 
it is same as the value observed for pristine phosphorene 
(0.17 mg). Note that, at higher fluorine concentration 
(x=0.833, 0.875), effective mass anisotropy is rotated by 
90°, such that m* in F-Y direction is almost two times 
larger than it’s value along F-X. For x=1.0, effective mass 
plot takes the shape of an ellipse, with major axis (i.e., 
higher m*) reverted to the F-X direction. In brief, due 
to fluorination m* calculated along the F-X direction de¬ 
creases, while it increases in the F-Y direction. Rotation 
of the effective mass anisotropy by 90° is observed only 
for x=0.833 and 0.875. 


IV. OXIDATION OF PHOSPHORENE 


Using cluster expansion method, at first we calculate 
the ECI values based on total energy calculation of 98 
structures using DFT method. A cross validation score 
of 20 meV indicates that the difference between fitted 
and actual (from DFT calculations) energy is negligi¬ 
ble. Next, using the calculated ECI values formation 
energies of 18650 phosphorene oxide structures (up to 32 
atoms per unit-cell) are predicted [see Fig. 3(a)] and the 
true and predicted “ground states” match very well. In 
Fig. 3(d), we illustrate the atomic arrangement and unit 
cell of the seven predicted “ground states” of phospho¬ 
rene oxide, namely P 12 O 2 , P 8 O 2 , P 12 O 4 , P 4 O 2 , PsOe, 
P 12 O 11 and P 4 O 4 for oxygen concentration ot x = 0.167, 
0.25, 0.33, 0.5, 0.75, 0.917 and 1.0, respectively. As 
shown in Fig. 3(b), the binding energy decreases with 
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increasing oxygen concentration. Note that, the oxygen 
adatoms are more strongly bound to monolayer phos- 
phorene at any given concentration, compared to fluo¬ 
rine adatoms [see Fig. 1(b) and Fig. 3(b)]. Moreover, 
contrary to the case of fluorination, it is noticed that ox¬ 
idation does not bring about significant changes in the 
P-P bond-lengths than compared to the pristine values. 
In fully oxygenated phosphorene, the longest P-P bond- 
length is found to be 2.38 A, which is 6.25 % more than 
the pristine P-P bond-length. However, significant orien¬ 
tational distortions are observed with bond-angles rang¬ 
ing up to 125° for P-P-P triplets for x= 0.75 and 0.917. 

A detailed study of electronic properties of the newly 
found “ground states” of phosphorene oxide reveals sev¬ 
eral interesting facts. Similar to the case of fluorination, 
first we rule out the possibility of any magnetic ground 
state of phosphorene oxide by doing spin-polarized elec¬ 
tronic structure calculation. Next, we find that with 
increasing oxygen concentration, bandgap initially in¬ 
creases from 0.9 eV in pure pristine phosphorene to 1.68 
eV at x=0.5, followed by a continuous decrease to 0.39 
eV in fully oxidized phosphorene [see Fig. 3 (c)]. The 
structure and the corresponding electronic band struc¬ 
ture of all the predicted “ground states” of phosphorene 
oxide are shown in Figs. 3 (d)-(e). As shown in the fig¬ 
ure, most of the “ground states” of phosphorene oxide are 
direct bandgap semiconductors (other than x=0.33 and 
x=0.75), unlike indirect bandgap ground states of fluori- 
nated phosphorene [see Fig. 1 (c)]. It is clear from the 
PDOS plots [shown in Fig. 3 (f)] that for x < 0.5, both 
the conduction and valence band states are dominated 
by orbitals belonging to the phosphorus atoms. How¬ 
ever, for X > 0.5 the valence band states are prevailed 
by the orbitals belonging to the oxygen atoms, while the 
conduction band states are still dominated by the phos¬ 
phorus orbitals. This is unlike fluorinated phosphorene, 
where both the valence and conduction band edges are 
dominated by the electronic orbitals of phosphorus at any 
given concentration [see Fig. 1 (f)]. 

The hole effective mass shows a very interesting be¬ 
havior after oxidation (see Fig. 4). In the F-X direction 
it decrease monotonically with increasing oxygen concen¬ 
tration; from 8.64 me in pristine phosphorene to 0.10 mg 
in fully oxidized phosphorene. On the other hand, in the 
F-Y direction m)( increases continuously from 0.16 mg in 
pristine phosphorene to 9.22 mg in fully oxygenated phos¬ 
phorene. Similarly, the electron effective mass, calculated 
along the F-X direction, decreases monotonically with in¬ 
creasing oxygen concentration; from 1.28 mg in pristine 
phosphorene to 0.08 Wg in fully oxidized phosphorene. 
On the contrary, in F-Y direction it increases from 0.17 
TOg in pristine phosphorene to 7.5 rUg at x = 0.75, be¬ 
yond which it decreases continuously and reaches a value 
of 0.4 mg of fully oxidized phosphorene. Interestingly, 
both hole and electron effective mass contours undergo 
a 90° rotation after oxidation, which takes place for x> 
0.75 and x> 0.25 for m)) and m*, respectively. 

Thus, the unique anisotropic electron and hole effec- 



FIG. 4. Polar plot of the hole and electron effective mass for 
various “ground states” of oxidized phosphorene. Note that, 
while both and m* mostly maintains a two-fold rotation 
symmetry, the anisotropy of the effective mass undergoes a 
90° rotation at higher oxygen concentration. 


tive mass of monolayer phosphorene can be modified in 
a controlled manner by adatom adsorption. This pro¬ 
vides desirable tunability in terms of device applications, 
because spatial conductance profile of phosphorene can 
be engineered to undergo a 90° rotation of the preferred 
conduction direction of both electron and hole. Compar¬ 
ing Fig. 2 and Fig. 4, we conclude that oxidation is more 
effective than fluorination for this purpose. 

Finally, let us discuss the choice of dangling vs. inter¬ 
stitial bridge position for the adsorption of oxygen atom. 
Since chemisorption in dangling location leads to minimal 
lattice distortion,it is expected to be the favorable po¬ 
sition for oxygen adsorption. Although this is found to be 
correct upto ^25% oxygen concentration, beyond which 
interstitial bridge position is predicted to be energetically 
more favorable.Thus, some of “ground states” found 
for higher oxygen concentration in this paper are possi¬ 
bly metastable states. Interestingly, it is found that the 
activation barrier for the surface oxygen atom (at dan¬ 
gling position) to move to an interstitial site (at bridge 
position) is ~0.7 eV.^^ Because of this, we believe that 
the dangling configuration could be kinetically stabilized 
even when they are not preferred thermodynamically. 


V. THERMAL STABILITY 

Before concluding, we confirm the stability of the pre¬ 
dicted “ground states” at higher temperature by doing 
Car-Parinello molecular dynamics^® simulations, as im¬ 
plemented in the quantum espresso package. Time step 
for MD simulation is taken to be 0.075 fs, with a total 
runtime of 1 ps and the fictitious electron mass (/r) is 
chosen to be 100 a.u to ensure that the electrons remain 
close to the ground state and good control of the con¬ 
served quantities is maintained. All the reported “ground 
states” are found to be stable well beyond room tern- 
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FIG. 5. The radial distribution function g{R) as a function 
of the radial distance R for T = 0 and T = 600K for (a) 
PFo .5 and (b) PO 0 . 5 . Snapshots of the structures are also 
illustrated in the inset. The broadening of the T = 0 sharp 
lines at finite temperature (600 K), due to thermal vibrations 
of the molecules is evident. 


perature. We test the stability in the range 400-700 K 
by checking whether the structures disintegrate due to 
thermally induced vibrations. We find that most of the 
“ground states” are stable upto 700 K, other than 
for X > 0.75 coverage, which are stable only upto 400-500 
K temperature range. 

Thermal stability for PF 0.5 and PO 0.5 is shown in 
terms of radial distribution function in Fig. 5(a) and 


(b), respectively. The sharp lines observed for 0 K are 
broadened at higher temperature due to thermally in¬ 
duced vibrations of the constituent atoms. We observe a 
peak broadening, measuring approximately ± 10 %, which 
is very similar to the values reported for graphene. 

VI. SUMMARY AND CONCLUSION 

To summarize, we systematically study chemical func¬ 
tionalization of monolayer phosphorene by surface ad¬ 
sorption of fluorine and oxygen. Using CE method to 
analyze all possible structures (upto 32 atoms per unit 
cell), we have found the ground states at different concen¬ 
trations of fluorine and oxygen, adsorbed on monolayer 
black phosphorene. It is observed that the bandgap"^^ 
as well as the effective mass of the charge carriers can 
be tuned by changing the concentration of adsorbates. 
In case of fluorination, most of the ground states are 
indirect bandgap semiconductors, while majority of the 
oxide “ground states” are direct bandgap semiconduc¬ 
tors. Other than few selective adsorbate concentrations 
(x=0.25 and 0.5 for fluorination and x=1.0 for oxidation), 
the bandgap of the newly predicted ground states are ei¬ 
ther comparable or higher in magnitude than that of pris¬ 
tine phosphorene. Further we find that the anisotropy of 
both the electron and hole conductivity of phosphorene 
can also be tuned via adatom adsorption, which can lead 
to the rotation of the preferred conduction direction by 
90°. From an application point of view, this phenomenon 
can be used to control and switch the easy direction of 
electronic and thermal transport in phosphorene based 
devices. 
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